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1. JNIRQaUCJIQN. 


Thla progress report for NASA Grant NGL 50-002-1 27 » Supplement No. 9 
"Multidisciplinary Research on the Application of Remote Sensing to Water 
Resources Problems,” at the University of Wisoonsin-Madison is for the pe- 
riod from August 1, 197B to July 31 » 1979. The Director of the Environ- 
mental Monitoring and Data Acquisition Group (now the Environmental Remote 
Sensing Center) , Professor Ralph W. Kiefer served as the Principal 
Investigator. The research was administc-.’ed through the Institute for En- 
vironmental Studies. 


Research on the application of remote sensing to problems of water 
resources was concentrated, beginning in 1976, on sediments and associated 
non-point-source pollutants in lakes. That work, except for one aspect, 
was concluded in this project year. Emphasis has shifted to further 
transfer of the technology of remote sensing to agencies which would use 
it, and to the refinement of equipment and programs for thermal scanning 
and the digital analysis of Images. The use of technology developed on 
the NASA remote sensing project by associated research projects has con- 
tinued this year to be an important contribution to many fields. 


In the area of instruction, with the approval as an official graduate 
program at the University of Wisconsin-Madison, of the Master of Science 
and Doctor of Philosophy in Environmental Monitoring last year, many stu- 
dents are pursuing these programs. The development of these programs and 
their growth to this stature have been directly dependent on the stimula- 
tion provided by this research project. 


As this project has completed its final year at full funding and en- 
ters the two years of reduced funding preparatory to its conclusion, ef- 
forts have been directed toward completing and documenting the development 
of techniques, assessing the future value of the techniques that have 
emerged, and preparing facilities and literature of value in tbs future 
application of this technology to new research problems. 


THE STRUCTURE OF THIS REPORT 


Section II. is the report of research projects funded wholly or in 
part by the NASA Grant. 

Section III. contains the reports of associated projects: those 

projects that depend on techniques developed on the NASA remote sensing 
grant, or that have made contributions of techniques to this research 
project. 


Section XV. has the Institutional Information about remote sensing 
research and Instruction at the University of Wisconsin requested by the 
Technical Officer, 

In Section V. publications emerging from the project, or closely re- 
lated to it are listed. 

The Appendix contains a paper that has depended upon basic techniques 
developed on the NASA remote snslng project. 


II. PROGRESS ON FUNDED PROJECTS 


/Lta9gCll?Ej.O. 


Professor Frank L. Scarpace 
James Verdin, Research Assistant 


A project Involving the use of multi-band Imagery to correct Landsat 
Imagery for atmospheric effects has been proposed for a number of years. 

In May of this year(1979) photographic Imagery was acquired over Green Bay 
Wisconsin coincident with a Landsat passover. The imagery was developed 
and found to be acceptable for the project. Since the Landsat Imagery was 
not available from EROS until Fall of 1979, however, no analysis of the 
data was possible during 1978-79. The student involved in the research 
has developed a method of overlaying a digitized photograph with the 
Landsat, but is leaving the University in August of 1979. It is antici- 
pated that another student will be involved in this project to complete 
the photographio/Landsat correlation during the 1980/81 academic year. 


Professor Frank L. Scarpace 
Linda Kalman, Research Assistant 


The goal of this portion of the research was to determine how 
photographic imagery and, ultimately, satellite imagery are affected by 
various types and concentrations of sediment in lake water. The specific 
goals of this sub-project during the year were to (1) finish the film 
calibration procedures and (2) start to analyze mission imagery acquired 
the previous two summers. The first of these goals was accomplished and 
the second task was started. 


Absolute calibration of film imagery was necessary to correct for at- 
mospheric problems and to use the imagery for radiometric calibration of 
Landsat imagery. Previous work had allowed us to calibrate the lens for 
geometrical lens falloff.* We were next interested in determining the ab- 
solute relationship between the energy falling on the film and the rela- 
tive exposure determined through normal sensltometrlc means. To this end, 
pictures of a calibrated lamp source were taken using several film and 
filter combinations. 


After development and scanning with our Optronic PI 700 scanning 
densitometer, average relative log exposure values were calculated for 
each image using existing computer programs. The average log exposures 
were calculated using training set programs which consisted of small sub- 
sets of pixels in the center of the lamp image. A calibration constant, 
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K| was calculated for each lamp Image which Is related to the relative ex- 
posure, and energy spectrum of the lamp, e(X). 


where Q is the p<'lid angle the lens subtends 


T is the exposure time 

F(l) is the transmittance of the filter used 
S(l) is the sensitivity of the film used 
is the "elA\cive log exposure 

e(X) is the energy spectrum of the calibration lamp 

Ideally this K should be independent of all parameters and should only 
vary a small amount from lens-to-lens. 


The experimental results seem to reflect this assumption with a de- 
termined K value for our system of 
K = 1.51 + .03 


There were a few data points as much as .1 larger or smaller than the 
above value. 


To check on the accuracy of the calibration constants, the absolute 
log exposures were calculated for images of reflectance panels. The imag- 
es were taken at the same time as incident flux radiometer measurements 
were made on the ground. 


The imagery used for analysis was photographs of reflectance panels 
taken in September 1978 using narrow band filters. The imagery was ac- 
quired at a flying height of approximately 700 m AMT on a clear to light 
haze day. While the imagery was acquired, incident flux measurements were 
performed on the ground using the digital radiometer with the cosine de- 
tector. The reflectance of the panels placed on the ground was measured 
the next day with the digital radiometer. Three shades of panels were 
used; white, dark gray, and light gray. The imagery was scanned at 50 
micrometers spot size. Training sets were established in the centers of 
the panel images and the program STATS was used to extract the average 
relative log exposure for the panel. The log exposures were corrected for 
lens falloff. The appropriate calibration factors found previously were 
applied to the corrected log exposures to yield absolute log exposures. 

The values were then converted to exposure, and divided by the product of 
w*t*AX to yield the actual intensity (ergs/ cm nm.sr.sec) incident on the 
lens surface. 


p 


E(X) incident on filter * 10^®®^^abs^*(^^^^^)erg8/eni‘‘"nra,8r.8eG 


Assuming no atmospheric effects, and a Lambertian type panel reflectance, 
R, the light incident on the panel is 

ft 

f 2 

r E(l) inc. on panel ■ E(X) inc. on filter * 1/R(X) ergs/era nratsr.sec. 


At this point, an average value of E incident on the panel was calculated 
for each wavelength, by averaging the results from all panels. Relatively 
small variation in calculated exposures for a given wavelength from one 
panel to the next was found. 

Since the panel was assumed to have a Lambertian surface, then the 
flux incident on the panel is merely 7 t<»E(X) incident on the panel. This 
value was computed for the three filters. 


f The flux incident on the detector radiometer (ergs/cm nm.sec.) was 

[ computed using the measured radiometer data and a calibration established 

earlier for the radiometer. The radiometer is equipped with a cosine de- 
tecto*^ (flat diffusing surface) for measuring incident flux from the 
I hemisp. ^ . a above the detector. The detector response was calibrated by 

f placing the flat surface of the detector flush against the surface of the 

’ standard lamp face. The standard lamp has an intensity calibration 

(ergs/cm nra.sr.sec.). The relative exposures calculated in this fashion 
t were in agreement with the measured fluxes by the radiometer within the 

expected accuracy of the film detector, 

I 

I 

\ The general conclusion of this part of the study is that the absolute ac- 

i curacy of an average of a number (N>3) of relative exposure determinations 

from film imagery is about 10%. This is consistent wich the log-detector 
character of the film. 


*Kalman, L. and F.L. Scarpace, ''Determination of Lens Falloff 
Through Digital Anlysis of Photographic Imagery," submitted to 
Photogrammetric Engineering and Remote Sensing, 
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Sail Simian SJiudlea. 


Professor Ralph W. Kiefer 
Donald R. Morris-Jonea 


Since the last progress report (which covered work completed up to 
about January of 1979) > a number of tasks have been completed by Donald R. 
Morris-Jones, a doctoral candidate in Environmental Monitoring supported 
on this project through May of 1979. An abstract of the dissertation to 
be completed by Mr. Morris-Jonea follows: 


ABSIRACI 


The Universal Soil Loss Equation (USLE) is a frequently used metho- 
dology for estimating soil erosion potential. A variety of types of 
geographic information (e.g. topographic slope, soil erodibillty, land 
use, crop type and soil conservation practice) are needed in order to 
quantify soil loss with the USLE. Topographic maps, soil surveys, field 
surveys and interviews with farmers are the traditional sources of this 
information. In this study, methods for gathering the land use/land cover 
Information required by the USLE were investigated with medium altitude, 
multi-date color and color infrared 70 mm positive transparencies and 
Landsat multispectral scanner data using human and computer-based inter- 
pretation techniques. In addition, the utility of computer geographic in- 
formation systems for storing, displaying and quantifying USLE data was 
also considered. Successful results, which compare favorably vfith tradi- 
tional field study methods, wero obtained within the test site watershed 
with airphoto data sources and human airphoto interpretation techniques. 
Computer-based interpretation techniques, which were applied to digital 
airphoto and Landsat multispectral scanner data, were not capable of iden- 
tifying soil conservation practices but were successful to varying degrees 
in gathering other types of desired land use/land cover information. For 
watershed and regional applications, computer geographic information sys- 
tems provided a number of significant advantages (e.g. synoptic coverage, 
retention of spatial detail, capability of simulating land use changes and 
their effects, identification of areas with severe erosion problems, etc.) 
over traditional, tabular methods of data storage, display and manipula- 
tion. 


The scope of the dissertation exceeds work completed under the NASA 
grant. 


At present, Mr. Morris-Jones is in the final stages of completing 
multidate interpretations of airphotos in order to derive cropping manage 
ment and support practice factor (i.e. C and P factors in the Universal 
Soil Loss Equation' values on a grid cell basis for the Pheasant Branch 
Watershed. These data will be stored in ;> detailed 'iomputer geographic 
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information systfsin of the watershed which has been developed to permit 
quantification of USLE values on a grid cell basis. 


Mr, Morris-Jones accepted a position as Research Engineer at the En- 
vironmental Research Institute of Michigan (ERIM) in the summer of 1979 
and is working within the Earth Resources Data Center. Landsat data proc- 
essing required for his dissertation will be conducted with ERIM facili- 
ties. Mr. Morris-Jones hopes to receive funding from contractors to de- 
velop and implement Landsat and other remote sensing based USLE informa- 
tion systems while at ERIM. Two proposals on this topic are presently un- 
der review. 


si£. a DIgit.al Imsa. Piosegsing 

Dr. Lawrence T. B'isher, Program Coordinator 


In our report for the 1977-78 project year we described the proposed 
acquisition of a McIDAS (Man-Computer Interactive Data Access System) col- 
or terminal. It is expected that this terminal •,'ill allow us to use mod- 
ern image- processing techniques which are now prohibitively expensive or 
slow. 


The McIDAS terminal will operate from a Harris/6 computer that the 
Meteorology Department is securing. It is hoped that new interactions be- 
tween information from meteorological satellites and that from Landsat 
will be developed as a fruitful area for research. 


A second component of our image-processing facility, the control com- 
puter for the Optronics P-1700 Scanning Densitometer, is also being re- 
vised. The Space Science and Engineering Center (SSEC) is building an In- 
tel microcomputer for control of the densitometer and tape drive to im- 
prove reliability and speed of the scanning process. 


The technicians of SSEC are now constructing the densitometer con- 
troller and the McIDAS terminal, with delivery expected in September 1979. 
Shakedown and programming will commence in the fall, with full operation 
predicted by December 1979. 


Conversion of programs for image analysis from the Univac 1110 format 
to that for the Harris/6 will be a formidable task, but completion of it 
should open many new areas of application of remote sensing for investiga- 
tion. 


j^cdirigafclgna iUi£ IHsrjal Saarmfir. 

Dp. Lowrenao T. c'ialier 


During tnla past year the data aequlaition syatera for our TI RC-ldA 
thermal acanner ayatem haa been changed from an analog recording system to 
an all digital recording ayatem. The Sangamo Saber III tape recorder naa 
been replaced by a micro-proceaaor (Intel 80/10) and a Pertee 1600 BPI 
tape deck. The acanner video output la tranaformed Into IP bit digital 
values (in order to record the full dynamic range of the poaaible signal), 
then recorded in ^IK byte records on the Pertee tape drive. 


The modifications were field tested on Lake Superior in the summer of 
1979 . As far as could be determined, the system is working up to expecta- 
tions. Future modification of the thermal scanner system will include a 
real-time video display of the imagery and the capability of some data 
analysis in the field. 


1 


J 
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III. ASSOCIATED PROJECTS 


fianapa.traftloa £tiQ.1sg,fr. EyaluiL&s, jai£ gpBQii&.unlslsa ficgriiBi.gr.fl]. 
AopIlaaBlana aJi Lapdaa.B fiata In Miaafloain 

Profeasor Frank L. Searpaoe 
Robert leadon, Research Assistant 


The objective of this study, beginning in Harch 1979 » was to demon- 
strate the use of Landsat data for the measurement of trophic status of 
I Inland lakes, in an operational setting in close cooperation with the user 

L agency. The Wisconsin Department of Natural Resources, as part of its 20b 

f Water Quality Program, must establish the trophic status of a large number 

of lakes. In this demonstration project, DNR scientists would operate the 
data e.'ctraotion and analysis programs developed at the University of Wis- 
consin, in collaboration with University research scientists. 

Original plans called for acquisition of Landsat computer-compatible 
tapes, extraction of data on the lakes in question, and analysis of the 
data— producing a list of lakes with trophic status indicated. Funds 
available for this demonstration, however, were adequate only for the ac= 
quisition of tapes and the commencement of data extraction. Analysis 
‘ would have to wait on further funding becoming available. 

Data have been extracted, and the project is ready for the next step 
f as soon as funds become available. 


F li&U aiudy 

I Prv/fessor Frank L. Scarpace 

\ Thomas Mace, Research Assistant 


The study of soybean crop yield using medium scale color and color 
infrared imagery to evaluate yield was completed. Methods were reported 
in our previous progress report (1977-78). Analysis was completed this 
year, and a final report prepared for Argonne Laboratory. A copy of the 
report is included as Appendix I to this report. 
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.Gaiai&firrAJLdffil filaaairisafrAflfl siL .Ofly-gr. 

Tncmas Maoe, KeaearoU Assistant 


Th3 Department of Forestry of the University of Wiseonsin, with the 
cooperation of the Environmental Remote Densing Center, oondueted this 
study to eompare forest cover type classification py digital analysis of 
small-scale (1:120,000) aerial Imagery to hand-drawn cover-type mapping 
from ground observation. 

The study area was Bear Island, Wisconsin, in the Apostle Islands Na- 
tional Lakeahore, Ground data were gathered by stratified random sam- 
pling. Aerial image classification was done by converting a color 
infrared aerial photo with density wedge to digital form, and classifying 
types by computer programs baaed on maximum likelihood with threshold. 

Results from the comparison of methods showed an average accuracy of 
62‘^ in the classification made by the computer-aided analysis of the 
small-scale aerial image. Greater accuracy could be expected if the orig- 
inal aerial data had been better adapted to this method. Imagery for op- 
timum use of this method of classification should have the following char- 
acteristics ; 

Scale 1:50,000 bo 1:100,000. 

The density wedge should be attached to and processed with the film 
image. 

Fall-off oliaracteristioa of the lens used for the original photo and 
filters used should be known in detail. No anti-vignetting fil- 
ters should be used. 

Images should be framed so that features to be classified are cen- 
tered. 

White reflectance panels should be included in the image when photo- 
graphed . 

If this method is used, and the suggested standards for imagery are 
met, greater accuracy can be expected, and maps of forest cover type can 
be made at substantially lower cost than by hand-drawing. This would 
greatly aid the National Park Service, which must map forest cover types 
over extensive areas in order to fulfill its goal of maintaining or recon- 
structing presettlement vegetation within National Parks and Monuments. 

The regular updating of these maps would also be accomplished more easily 
and at less cost. 
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firaan £r.Q.1agA, 




Professor Frank L. Scarpace 
Bruce Quirk, Research Assistant 


This project Involves development of techniques for mapping lane cov- 
er from Landsat and digitized aerial photography for use In a hydrological 
runoff model. It Is anticipated that the resultant data will be used by 
the Wisconsin Department of Natural Resources In their 208 Water Quality 
programs as well as by the United States Geolrglcal Survey for Input to 
their streamflow-water quality models. 


Thus far, Landsat and aerial Imagery have been acquired for the test 
site. The aerial Imagery has been converted to digital form and classifi- 
cation of test sub-watersheds has begun. The Input of the land cover data 
to the sediment run-off models will be started during the next year. 


This research Is funded by the University of Wlsconsln-Madlson Sea 
Grant program. 


fjioaganfr, Bcangh watershed Soli Study. 

Professor Ralph W. Kiefer 
Akin Fapohunda, Graduate Student 


The general perception in literature in a sense recognizes some limi- 
tations as to the degree of detail and accuracy that can be obtained using 
remote sensing. Thus It is agreed, for instance, that soil mapping at the 
association level can adequately be obtained even from spaceborne platform 
sensors such as Landsat. What is generally overlooked is that while this 
can be achieved in some places, it becomes increasingly difficult and per- 
haps impossible in others. For example, in terrains with relatively young 
geological and pedalogical age such as wetlands, the characteristics on 
which photographic processes are based are not fully expressed. Attempt- 
ing the application of remote sensing in these areas is at best only pos- 
sible with such a detailed level of field work to gather ground truth that 
whatever advantage results is eventually minimized. 


The aim of this research is to demonstrate that the level of field 
work required for soil mapping using aerial photography varies spatially 
depending on terrain geology and complexity of soil association systems. 
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Three scenes, all in the Pheasant Branch Watershed in Dane County, 
Wisconsin, were selected for the purpose of this study. The interpreta- 
tion to delineate recognizable soil bodies was done on the color additive 
viewer and zoom stereoscopes. Overlays of delineations obtained in the 
area covered by each photograph were prepared for field use. 


Field work was carried out for four weeks during summer 1979* The 
results indicate significant variations in the amount of information con- 
tained in each photograph and consequently the amount of time required for 
field work to produce accurate soil classification using remote sensing 
techniques (in this case, color and color infra-red photography). 


A detailed report of this work is being prepared for publication by 
May, 1980. 


^apjJLpg sL Vegetation iii sng.bgygan idaTjaii 

Professor Ralph W. Kiefer 
Professor Frank L. Scarpace 
Bruce Quirk, Research Assistant 


Since 1975 the remote sensing group has been studying methods for 
mapping vegetation in the Sheboygan Marshy a wetland in central Wisconsin. 
Small-scale imagery (1:120,000) has been used to determine how useful it 
is for classification of vegetation in the wetland and how vsll the bound- 
aries of the wetland can be delineated. 

This study has been completed. A full report will be included in the 
next progress report. 


Data g Ql la ff tl gn . aM Analysis Methods Used M Monitor Impacts 
Qz?r. Ujaa la a severely Disturbed Wetland 

Professor Ralph W. Kiefer 
Sarah L. Wynn, Research Assistant 


The purpose of this study has been 1) to document wetland species and 
community change from 1975 to 1977 from the air and from the ground, and 2) to 
evaluate and compare the efficiency, sensitivity, and reliability of a variety 
of traditional and new data sampling and analysis techniques used to document 
this change. This study was part of a larger, ongoing Environmental Protec- 
tion Agency study to monitor the impacts of the construction and operation of 
a coal- fired power plant located in a large wetland near Portage, Wisconsin. 




The nine methods chosen for this study were evaluated on a basis of time, 
cost, sensitivity, and reliability and also on the bases of: whether they 

monitor vegetation change at a species or community level,- whether they moni- 
tor community change in terms of change in area and/or location; and whether 
they provide information about community trends. 


Four of the nine methods (diversity index, subjective classification, as- 
sociation analysis, and structure analysis) used ground sampling data. Two 
methods (airphoto monitoring and disturbance mapping) used airphoto data only, 
while airphoto grid analysis and airphoto interpreted vegetation mapping used 
airphoto and ground sampling data. The ninth method (computer-assisted map- 
ping) used only airphoto data but relied heavily on the analyst’s knowledge of 
the area. 


A Ph.D thesis based on this work will be completed by December 1979. 


COLUMBIA (lEIiEMmifiL STAIIQM. iffilLAlffiS. iSfiETAIIOR ,QB.QB£ 

Barbara Bedford 


The Wetlands Plant Ecology Group of the Columbia Generating Station Im- 
pact Study has studied the effects of leakage from the station's cooling lake 
on wetland vegetation. The study included four phases: field inventory and 

classification of plant communities, vegetation monitoring, field and labora- 
tory experiments to test hypotheses regarding the mechanisms controlling popu- 
lation changes, and assessment of field and theoretical approaches. 


Results showed that changes in water levels and water temperatures caused 
by seepage from the cooling lake led to significant changes in wetland plant 
populations and communities within one year after the station began operation. 
In general, dominant perennial rhizomatous species of the genus Carex de- 
creased in density and distribution, hydrophytic species such as Tvoha 
latifolia increased, and annual species such as Bidens oernua and Pilea 
Dumila . which had been insignificant or absent before disturbance, increased 
markedly. 


Results of community-level measures of wetland response to leakage from 
the cooling lake showed significant changes occurring in the structure of the 
plant communities. Major shifts in dominance and diversity patterns and a 
continuing trend of decreasing vegetative cover occurred. Open water and ex- 
posed mudflats replaced the previously closed and densely vegetated perennial 
plant communities over an increasing portion of the study area. Annuals 
colonized some of the habitat opened by the removal of perennial species, but 
large areas remained unvegetated. The plant commenced operation in 1975. By 
1977 1 19^ of the quadrats sampled in the area of major impact had no rooted 
vegetation. Another contained only annual vegetation. 
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This study indicated that prediction of the impact of such a disturbance 
on wetland vegetation will require a knowledge of the key characteristics of 
the life cycle of the plants involved. 


Remote sensing of the wetland vegetation has been used by this group 
through analysis of medium-altitude aerial imagery acquired by the Environmen- 
tal Remote Sensing Center using the airplane of the Wisconsin Department of 
Natural Resources. Analysis of photos has allowed the group to; 1) observe 
changes through time in spatial patterns of the vegetation, 2) map patterns of 
change in vegetation, 3) delineate the boundaries of areas of open water, ‘i) 
determine areas of high leakage from the cooling lake, 5) compare seasonal 
alterations in vegetation pat'^i;;* s caused by warm water leakage with control 
site seasonal variations, 6! : t jte areas of high temperature water during 
spring, fall, and winter. 


Acquisition of data on the ground in this wetland is hazardous or impos- 
sible at certain times of the year due to the peculiar combination of water 
levels, floating aats of vegetation, and dense growth. 


The study is funded by the U.S. Environmental Protection Agency. 


IV. SUMMARY OF PROJECT ACTIVITIES 


PROJECTS COITDUCTED "./ITH NASA FUllDS 
PROJECT DATES 

ATtlOSPHERIC CORRECTIONS 

Simulationeous low-altitude aerial photographs and landsat imagery 1976-79 
have been acquired. Relations of these are analyzed tc determine 
a means for correcting satellite 

SEDIMENT PLUME STUDY 

Low and medium-altitude photographs of sediment plumes in an inland 1976—79 
lake were acquired in two successive years during runoff events. 

Progress is being made in calibrating the images to determine if 
sediment concentrations can be measured by remote sensing of this 
kind. 


SOIL EROSION S TUDIES 

Medium altitude aerial photographs were analyzed to secure landcover 1977-79 
for use with other factors in the Universal Soil Loss Equation to 
determine erosion due to runoff. Measurements by remote sending 
are compared to those made by ground measurement. 

SCANNING DENSITOMETER CONTROLLER AND McIDAS COLOR TERMINAL 

A new computer and tape drive for use with the P-1700 Scanning Densi- 1978-79 
tometer and a color McIDAS display and analysis terminal are being 
acquired to allow quicker and more precise image processing for both 
aerial photographs and Landsat images. 


MODIFICATIONS TO THE THERMAL SCANNER 

The Texas Instruments RS-18A Thermal Scanner is being modified to 
record signals in digital rather than analog form. 


ACTION/DECISION 


Data acquired and prepara- 
tions are being made for 
analysis 


Film being calibrated. 
Sediment gradients are not 
readily measurable by this 
means . 


Remote sensing provides a 
practical means to deter- 
mine some factors for the 
USLE, 


Equipment is now being 
constructed 


Scanner is in operation. 
The reliability of 
operation seems much 
improved . 


ASSOCIATED PROJECTS 


These projects share technology with the NASA project, but are funded from other 


PROJECT DATES AGENCY FUNDING 

GREEN BAY PROJECT 

Development of techniques for mapping land cover 1978-81 Sea Grtaiit 
from Landsat imagery and digitized aerial photo- 
graphs for use in a hydrologic runoff model. 

Results may be used by the Wisconsin DNR in 
their 208 Water Quality Program, and by the USGS 
in streamf low-water quality models. 

PHEASANT BRANCH WATERSHED SOIL STUDY 

This study seeks to demonstrate that the level 1978-79 Self-funded 

of field work required for soil mapping using by student, 

aerial photographs varies spatially depending 
on terrain geology and complexity of soil 
association systems. 


COMPUTER-AIDED CLASSIFICATION OF FOREST COVER 
TYPES FROM SMALL SCALE AERIAL PHOTOGRAPHY 
High-altitude color infrared photographs of 
forests were analyzed to determine how 
accurately and at what costs forest cover 
types could be determined. 


1978-79 University of 
Wisconsin Dept, 
of Forestry. 


DATA COLLECTION AND ANALYSIS IIETHODS USED TO 
MONITOR BIPACTS IN A WETLAND 
Remote sensing by means of medium-altitude, 
color infrared aerial photographs is compared 
with ground surveys for monitoring changes 
to vegetation in a wetland. 


1975-79 U.S. Environ- 
mental Protection 
Agency 


sources 

ACTION/DECISIQN 


Data acquired analysis 
and model study under 
way. 


Significant variation 
found in the amount of in- 
formation discernible from 
different photos and con- 
sequent differences in the 
field work required for 
accurate mapping. 


Accuracy of 62%, with 
many improvements sugges- 
ed. Cost not calculated. 


Remote sensing is recom- 
mended as a tool for 
wetland monitoring 


ASSOCIATED PROJECTS 


These projects share technology with the NASA project, but are funded from other sources 


I 

I 


PROJECT DATES AGENCY FUNDING 

COLUMBIA GENERATING STATION 
WETLANDS VEGETATION STUDY 

Color infrared aerial photographs of a wetland 1975-78 U.S. Environ- 
ed jacent to a power plant cooling lake were used mental Protection 

to document vegetation community changes in Agency, 

the wetland. 


ACTION/DECISION 


Remote sensing is 
a valuable technique 
providing a time- 
series record, and 
access to difficult 
areas for study. 


MAPPING OF VEGETATION TYPES IN SHEBOYGAN MARSH 

From high altitude (1:120,000 scale) color 1978-79 University of 

infrared Imagery, maps of the vegetation Wisconsin 

types in a wetland are being prepared and 
compared with field data. 


Full report and 
correlations will be 
included in the 
1979-80 progress 
report . 


DEMONSTRATION PROJECT TO EVALUATE THE 
OPPORTUNITIES FOR OPERATIONAL APPLICATIONS 
OT LANDSAT DATA IN WISCONSIN . 

In cooperation with the Wisconsin Department 1979-80 NASA-ERRSAC 

of Natural Resources, this study seeks to 

classify lakes in Wisconsin by trophic 

status using Landsat. Classification is 

done by DNR personnel under supervision 

by University research personnel. 


Funds for only a start 
in data extraction. 
Will seek more funds. 


PROJECTS LISTED BT COOPERATING AGENCI 


Es^nal 


U. S. Geological Survey 
SEDIMENT PLUME STUDY 


U. S, Environmental Protection Agency 

DATA COLLECTION AND ANALYSIS METHODS USED TO MONITOR IMPACTS IN A WETLAND 
COLUMBIA GENERATING STATION WETLANDS VEGETATION STUDY 


NASA - Eastern Regional Remote Sensing Applications Center 

DEMONSTRATION PROJECT TO EVALUATE THE OPPORTUNITIES FOR OPERATIONAL 
APPLICATIONS OF LANDSAT DATA IN WISCONSIN 


State 


Wisconsin Department of Natural Resources (DNR) 

SEDIMENT PLUME STUDY 

ATMOSPHERIC CORRECTIONS 

SOIL EROSION STUDIES 

MODIFICATIONS TO THE THERMAL SCANNER 

MAPPING OF VEGETATION TYPES IN THE SHEBOYGAN MARSH 

DEMONSTRATION PROJECT TO EVALUATE THE OPPORTUNITIES FOR OPERATIONAL APPLICATIONS 
OF LANDSAT DATA IN WISCONSIN 


Wisconsin Geological Survey 

SEDIMENT PLUME STUDY 
SOIL EROSION STUDIES 
GREEN BAY PROJECT 

PHEASANT BRANCH WATERSHED SOIL STUDY 


Rgglgnal 


Wisconsin Power and Light Company 

DATA COLLECTION AND ANALYSIS METHODS USED TO MONITOR IMPACTS IN A WETLAND 
COLUMBIA GENERATING STATION WETLANDS VEGETATION STUDY 




City of Middleton 
SEDIMENT PLUME STUDY 


PROJECTS LISTED BY SOURCE OF IMAGERY 


Landsat Imagery 
ATMOSPHERIC CORRECTIONS 

SCANNING DENSITOMETER AND McIDAS COLOR TERMINAL 
GREEN BAY PROJECT 

DEMONSTRATION PROJECT TO EVALUATE THE OPPORTUNITIES FOR OPERATIONAL APPLICATIONS 
OF LANDSAT DATA IN WISCONSIN 


High Altitude Imagery 
GREEN BAY PROJECT 

COMPUTER-AIDED CLASSIFICATION OF FOREST COVER TYPES FROM SM.^'.L-SCALE 
AERIAL PHOTOGRAPHY 

MAPPING VEGETATION TYPES IN SHEBOYGAN MARSH 


Medium Altitude Imagery 

SEDIMENT PLUME STUDY 
SOIL EROSION STUDIES 
PHEASANT BRANCH WATERSHED SOIL STUDY 

DATA COLLECTION AND ANALYSIS METHODS USED TO MONITOR IMPACTS IN A WETLAND 
COLUMBIA GENERATING STATION WETLANDS VEGETATION STUDY 


/•; 
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LOW Altitude Imagery 

SEDIMENT PLUME STUDI 
ATMOSPHERIC CORRECTIONS 


OTHER FUNDING RECEIVED 


These amounts, related to the capabilities generated by this grant, 
are not auditable, since they represent portions of larger grants based on 
our estimate of the cost of that part related to remote sensing. 


.Esilsral fiRgaslg.8 

U. S. Environmental Protection Agency $35,000 

NASA - Eastern Regional Remote Sensing Applications Center $9,db0 


S.fca,fc& Agsnslsa 


University of Wisconsin-Madisen 
Wisconsin Department of Natural Resources 
Sea Grant Program 


$50,000 

$ 2,000 

$39,000 
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EDUCATIONAL ACTIVITIES 


u 


Professor Scarpace participated in a series of discussions with the Wis- 
consin Department of Natural Resources on plans for a program of classifi- 
cation of Wisconsin lakes using Landsat, under funding from NASA, ERRSAC. 

Professor Kiefer continued his work as a member of the State of Wisconsin 
Interagency Wetlands Mapping Advisory Committee. Plana for remote sensing 
activities for all State agencies concerned with wetlands mapping are co- 
ordinated by this committee. 


The following courses in Remrte Sensing were offered on campus, with en- 
rollments as shown: 


.Courgg, N-UPibgc 

.C.pjadl.t.a. 

CQUt:se Name 

Enrollment 

CEE/IES-552 

3 

Remote Sensing of the 
Environment 

42 

CEE 55 ^ 

3 

P'undamentals of Remote 
Sensing 

14 

CEE 555 

3 

Airphoto Interpretation 

46 

IES/CEE/LA-556 

3 

Remote Sensing Image 
Interpretation 

12 

CEE 351 

3 

Photogrammetry for Non-Engineers 

14 

CEE 356 

3 

Photogramraetry 

39 

lES/CEE 920 

3 

Environmental Monitoring Seminar 

15 


Numerous courses dealing with advanced photogramraetry, surveying, and car- 
tography are also offered. Master of Science and Doctor of Philosophy de- 
grees with an emphasis on Remote Sensing can be earned in either the De- 
partment of Civil and ]jJ;j/ironmental Engineering or in the 
Intfardisciplinary Environmental Monitoring Program. 
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AGENCY CONTACTC 


( 


F.ggsr.al 


Environmental Protection Agency; Dr. Gary Glass, Research Chemist, 

National Environmental Research 
Laboratory, Duluth, Minnesota 

Tennessee Valley Authority: Allan Vosa,RewoLo Sensing and Photogramme try 

Section, 

Chattanooga, Tennessee 

NASA - Eastern Regional Remote Sensing Applications Center: 

Betsy Middleton 
902.1, Building 22, Room G56 
Goddard Space Flight Center 
Groenbelt, Maryland 20771 

U. S. Geological Survey: Greg Allord, Cartographer 

1815 University Ave. 

Madison, WI 53706 


S.ka.ktt 


Wisconsin Department of Natural Resources; GEF 2 

Madison, WI 53702 

John M. Cain, Chief, Water Quality Planning Section, 

John G. Konrad, Chief, Special Studies Section, 

Jerome R. McKersie, Chief, Water Quality Evaluation Section 
Francis Schraufnagel , Chief, Water Quality Section. 

Wisconsin Geological Survey: 1815 University Ave. Madison, WI 53706 
Arthur Ziegler, State Cartographer 
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TKAFFIC TO THE ENVIRONMENTAL REMOTE SENDING CENTER 


Requests for Landsat Imagery 

Requests for aerial photographic imagery 

Contacts by citizens seeking infomation on remote sensing 

Inquiries from State and local agencies 

Consultations with other research or instructional staff for 
consideration of applications of remote sensing 


GRADUATES OF THIS PROGRAM 


Steve Fratoni, Natural Resource Data Management Study, 
New England Water Basic Commission, 
Boston, Mass. 

James Verdin, Water Resources Program 
Colorado State University 
Fort Collins, Colo, 

Donna Stetz, EROS Data Center 

Applied Science Division 
Sioux Falls, S.D. 
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PUBLICATIONS 


Members of the research team wrote or contributed to the following publi- 
cations during this project year. Some of these were listed in the previous 
progress report | but are included here because they came out during the porioa 
covered by this report. 


Allord, G,, and F» L. SoarpacOf "Improving Stream Flow Estimates Thru the Use 
of Landsat," Proceedings , AWRA Satellite Hydrology, June 1979. 

Eahart, J, , E, Joeres, and J, Hoopes, "Optimization of Waste Removal for Wide, 
Shallow Rivers," J.9,ur.n,al Sit £nYlR0nBgn.tal EnKln,fi,eriflJS BiYlaiQUi American 
Society of Civil Engineers, 197U, 

Fisher, L. T., F. L. Scarpace, and R. Thompson, "Multidate Extraction Proce- 
dures for Landsat Lake Quality Monitoring," f.hO.tOffiramBetrlO Engineer.img aM 
H, gag to Vol. 45, No. 5, May 1979. 

Green, T,, J. Hoopes, F. Scarpace, J. Villemonte, and R. Madding, "Measuring 
Effluent Plumes," Proceedings of the ASCE Conference on Hydraulic Research, 
August 197H. 

Green, T,, and R, Terrell, "Airborne Radiometer Measurements of the Keweenaw 
Current." Journal of GeoDhvsloal Research. 83, pp. 419-426, 197b. 

Kalman, L., and F. L. Scarpace, "Determining Lens Fall-off Through Digital 
Analysis of Photographic Imagery," Proceedings of the 45th Annual Meeting of 
the American Society of Photogramme try, March 1979. 

Llllesand, Thomas M., and Ralph W. Kiefer, "Remote Sensing and Image Inter- 
pretation", 1979, John Wiley and Sons, Inc,, New York, N.Y, 612 pages, 

Melnert, D., D. Malone, A, Voss, and F, L. Scarpace, "Trophic Classification 
of Tennessee Valley Area Reservoirs Derived From Landsat Multispectral Scanner 
Data," Tennessee Valley Autl r*lty Technical Report, November 1976, 

Morgan, Ken M. , Gerhard B. Lee, Ralph W. Kiefer, T. C. Daniel, Gary Bubenzer, 
and John T, Murdock, "Prediction of soil loss on cropland with remote sens- 
ing," ianKnal St gall 3ml Misn Ssns^ny^ilm, Nov, -Deo, 1976, V0l,33, Number 
6, pp, 291-293. 

Morgan, Ken M., Donald R. Morris-Jone.’, Gerhard P„ Lee, and Ralph W, 

Kiefer, "Cropping Management Using Color and Colo” Infrared Aerial 
Photographs," Photogrammetrlo Engineering and Remote Sensing . Vol, 45, No, 6, 
June 1979, pp. 769-774. 

Morris-Jones, D. R, , Ken M. Morgan, R, W, Kiefer, and F, L, Scarpace, "Remote 
Sensing as a Source of Land Cover Infoni'atlon Utilized in the Universal Soil 
Loss Equation," Thirteenth International Symposium on Remote Sensing of Envi- 
ronment, April 1979. 


A 


Scarpace, P, L., T. Green, and R. Madding, "A Thermal Scanning Study of 
Coastal Upwelling in Lake Superior," Proceedings of the 45th Annual Meeting of 
the American Society of Photogramraetry, March 1979. 

Scarpace, P. L., K. Holmquist, and A. Voss, "An Investigation of Atmospheric 
Corrections to landsat Data," P_rQceedings of the AWRA Satellite Hydrology Con- 
ference, June 1979. 

Scarpace, F. , "Densitometry on Multi-emulsion Imagery," Photogrammetrio Engl- 
neering and Remote Sensin g. Vol.44, No. 10, October 1978. 

Scarpace, F. L., and G. A. Friedrichs, "A Method of Determining Spectral Ana- 
lytical Dye Densities," EMing.s.nini?: ^ JiaiiKitg. &snslnR, 

Vol.44, No. 10, October 1978, pp. 1293-1301 . 

Scarpace, F. L. ,D. Meinert, D. Malone, and A. Voss, "Trophic Classification 
of Tennessee Valley Area Reservoirs," Proceedings of the 46th Annual Meeting 
of t..e American Society of Photogramraetry, March 1979. 

Scarpace, F. L., T. Green, and R. Madd'g, "A Thermal Scanning Study of 
Coastal Upwelling in Lake Superior," Proceedings of the 45th Annual Meeting of 
the American Society of Photogramraetry, March 1979. 

Scarpace, F. L., and T. Green, "The Spatial Variability of Coastal Surface 
Water Temperature During Upwelling," Journal of Physical Ooeanography. , Vol.9, 
No .3, May 1979. 

Scarpace, F. L,, K. Holmquist,, and L. T. Fisher, "Landsat Analysis of Lake 
Quality." Photogrammetrio Engineering and Remote Sensing . Vol.45,No.5, May 
1979. 
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REPORT TO ARGONNE NATIONAL LABORATORY: 


MIcrodcnsi tometrlc Analysis of Aerial Photographic Imagery 
for Yield Prediction In SO^ and Acid Rain-Damaged Soybeans 


' by 

Thomas H. Mace 
Frank L. Scarpace 


February 20, 1980 


iMAOfiUY ACQUISITION 
A(»rin I Ph otogrophy 


Aerlol photography was taken over the stuJy area on 22 August, August, 
and 6 September 1978. Three 70 mm llasselblod cameras were vertically mounted 
In a Cessna 180 aircraft. All cameras were equipped with 80 mm focal length 
lenses, and each mission consisted of passes at 610, 1220, and I830 meters 
above ground level. Camera 1 was loaded with Kodak 2Vi8 (color reversal) 
film and equipped with a haze filter; camera 2 was loaded with Kodak 2^3 
(color Infrared reversal) film and equipped with a yellow {H\2) filter; and 
camera 3 was loaded with Kodak 2l)03 (tri-x, panchromatic negative) film and 
no filter. All film had been stored frozen and was shipped to Precision Photo 
Lab Inc., Dayton, Ohio (via air courier) for sensitometry and processing. 

Upon receipt of the processed film, we discovered that the film from the 
22 August mission had been improperly exposed and v/as of marginal value. 
Apparently, the magazines holding the color and color infrared film had been 
switched, resulting in the color film being exposed through a yellow filter, 
and the color infrared film having been exposed through a haze filter. Sim- 
ilarly, processed tri-x was substituted for unexposed tri-x, making the 
imagery from camera 3 unusable as well. All other imagery appeared to be 
exposed properly and to be processed correctly. 

Scanning 

Table 1 lists the imagery scanned by: frame number, date of photography, 
film type, and nominal scale. All images were scanned on an Optronics PI700 
rotating-drum-type scanning microdensitometer at an aperature/rastcr setting 
of 5*^ pm. The data for each image, and its corresponding film wedge, were 
written on magnetic tape. 

\ 


r/ 


TABLE ]. 


IMAGERY SCANNED - ALL SCANNED AT 50u SPOT SIKE 


FRAME NUMBER 

DATE 

FILM 

NOMINAL SCALE 

62-/i57-25 

29 Aug 78 

Color 

1:7620 

62-ii57“3 

29 Aug 78 

Color 

1:301)80 

62-ii58-20 

6 Sep 78 

Color 

1:7620 

62-1)58-/1 

6 Sep 78 

Color 

1:301)80 

62-1)56-25 

29 Aug 78 

CIR 

1:7620 

62-1)56-3 

29 Aug 78 

CIR 

1 : 301)80 

62-1)55-21 

6 Sep 78 

CIR 

1 : 762.0 

62-1)55-1< 

6 Sep 78 

CIR 

1 : 301)80 

62-1)59-3 

29 Aug 78 

TRI-X 

1:301)80 

62-1)59-25 

29 Aug 78 

TRI-X 

1 :7620 

62-1)60-3 

6 Aug 78 

TRI-X 

1 : 301)80 

62-1)60-9 

6 Aug 78 

TRI-X 

1:7620 

*60-782-2 

22 Aug 78 

CIR 

1 : 301)80 

*60-782-22 

22 Aug 78 

CIR 

1 : 7620 

*60-783-1) 

22 Aug 78 

Color 

1 : 301)80 

*60-783-2!) 

22 Aug 78 

Color 

1 :7620 


The multi-emulsion films (color and color infrared) were scanned three 
times over the same area. One of three narrow band (10 nm bandwidth) 
filters was placed in the optical path, such that the three scans recorded data 
in the blue (.^5 pm), green (.55 pm), and red (.65 pm) regions of the spectrum 
These data are proportional to the integral dye densities of the three film 
layers (at the three measurement wavelengths) and are related to exposures 
from the field by the sens! toinetr ic film wedge and the characteristics of 
the film dyes (Scarpace, 1,978.). 

The single-emulsion images (tri-x), and their corresponding film wedges, 
were scanned once (through the 550 nm filter). The density/exposure rela- 
tionship, in this case, is expressed by the film v/edge alone, as no dye 
Interactions are present. 


Thus, the result of scimniny was that tlie inieges were transform'd 
Into diflitol arrays (1200 x 1200 x 3 fot* color and color Infrared, and 
1200 X 1200 X 1 for tri-x) of density measurements every 50 )im spatial 
locution on the film. These 50 pm measurements are coninonly called 
"picture cle.nents" (or "pixels") and consist of integers from 0 to 255 
corresponding to densities from 30 to OD, respectively. 


IMAGERY ANALYSIS 

Imagery Selection and Correction 

The images from the 29 August mission were selected for intensive 
analysis. The multidate aspect of the study was abandoned due to the 
poor condition of the 22 August imagery and the lack of adequate reference 
points in the 6 September Imagery.^ It was decided, therefore, to analyze 
the data v/ith respect to sensor (film type) and scale. 

Correction of the data sets for dye interactions and lens falloff were, 
performed according to the procedures outlined in Scarpace (1978) and 
Kalman and Scarpace (1979). The values stored for each pixel were, therefore, 
relative log exposures which had been scaled from 0 to 255- 
Data Extraction 

Line printer character maps of the corrected data arrays were made, and 
the locations of the field plots were outlined. The data points within the 
outlines were extracted and a mean vector for each plot was calculated. 

^The dispersal pipes were used for photogrammetr I c control and In locating 
the plots. These pipes had been removed by 6 September, and It was felt 
that we could not locate the plots as accurately on that date. Thus, a 
multidate evaluation might not be strictly valid. 


y V 
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These ftKians v;ere passed to another program whiclj transformed the scaled mean 
relative log exposures to relative exposures (unsealed, real numbers). These 
mean exposures were subsequently used for statistical analysis, 

Additionally, tv/o methods were used to achieve lumped means as well as 
individual plot means for each treatment. The first method was to simply 
calculate means by plot (le, ILC-l, H.C-2, UC-3, etc,) The second method 
involved merging the data from the four common treatment typos and then cal- 
culating the means. Further, Inspection of the large scale Imagery (1:7620) 
revealed that every other data row corresponded to a plant row. Therefore, 
pixels whose values Indicated they were non-plant pixels (le, shadow and soil) 
were removed from the merged plots before calculating the means. It should 
be noted, however, this could only be done for the 1:7620 scale imagery. 
Regression Analysis 

Tables 2 and 3 summarize the results of the regressions performed 
using the HINITAB program package on the University of Wisconsin Unlvac 1110 
computer. The terms "ALL,” "SO^," and "A.R. S SO^" refer to the two separate 
experiments taken together, then Individually. In each case, corresponding 
control plots were included. The expressions in parentheses indicate which 
spectral bands were used as predictors for yield. The term "(^,5,6)" was 
used when the three predictors were: 1) blue, green, and red bands for color 

film, ?) green, red, and Infrared bands for color infrared film, or 3) a 
single broad band over blue, green, and red for the black and white panchro- 
matic film. Method A refers to no merging or masking of plots, and method B 

refers to plots which had been merged and masked (1:7620 only). The values 
2 

of R ranged from 3*U to 91.6%, and the values of S ranged from 88.8 kg/ha 
to 2^7*0 kg/ha. 

The results were sometimes conflicting, but, generally. In accordance 


with whnt one would expect. Method B wa& ur.uolly an idiproveinetu over method A. 
This may be partially the result of more point c being used for each mean cal- 
culation and partially tne result of fewer points being used for the linear 
regression. Except for large scale acid rain experiment plots, 3 predictors 
were better than 2 ., and the individual band predictors followed in the 
order; red, green. Infrared, and blue. 

Soma interesting differences occur when one compares the regressions 
from the SO2 plots to the regressions from the acid rain plots. There was 
a sharp increase in tlie usefulness of the infrared data vdien the acid rain 
plots were analyzed. This was particularly true for the low altitude (large 
scale) data. This would indicate that acid rain changes tlie. infrared reflec- 
tance of the leaves in a way which Is more related to yield than does SO2. 

However, one must consider this phenomenon more an observation than 
a conclusion. We suspect that there may be interactions occurring between 
the variables v/hich are not completely represented by our data. For example, 
the general improvement (in all but the infrared regressions) in R for the 
acid rain plots with an increase in flight altitude is perplexing. 
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R^% 


Film Type 

COLOR 

COLOR IR 

TRI-X 

COLOR 

COLOR IR 

Frame: 

457-25 

456-25 

459-25 

457-3 

456-3 

Plots (Bands) 







ALL(4,5,6) 

68.9 

81.0 

54.8 

65.0 

56.1 59.3 

58.9 

73.7 

69.2 

SG.4 

(G,R) 

53.4 

63.1 

52.9 

62.fi 

- - 

58.8 

73.2 

62.6 

74.8 

(B) 

19.3 

8.1 

- 

- 

- 

36.9 

38.2 

- 

- 

(B) 

43.5 

57.0 

48.9 

47.2 

- 

53.1 

62.1 

56.5 

71.0 

(R) 

- 

- 

52.9 

62.2 

- 

53.9 

70.7 

60.4 

71.3 

(IR) 


- 

14.5 

48.0 


- 

- 

55.8 

63.7 

502(4,5,6) 

70.5 

91.6 

53.3 

81,3 

71.4 86.8 

54.3 

72.1 

67.5 

91.6 

(G,R) 

53.7 

64.7 

52.7 

74.3 

- - 

54.3 

72.1 

59.4 

73.6 

(B) 

19.2 

3.1 

- 

- 

- ~ 

37.4 

52.6 

- 

- 

(G) 

44.6 

61.1 

48.8 

71.4 


49.4 

65.2 

53.1 

73.0 

(R) 

53.5 

64.4 

52.6 

69.3 

- 

49.6 

63.1 

57.S 

68.4 

(IR) 

- 

- 

11.5 

38.5 


- 

- 

54.2 

63.5 

A.R.(4,5,6) 

74.2 

78.8 

77.4 

89.3 

50.5 32.5 

75.5 

75.2 

80.6 

£0.4 

a (G,R) 

59.3 

60.3 

60.8 

56.2 

- - 

74.8 

74.4 

£0.5 

£0.2 

SO 2 (B) 

12.5 

16.5 

- 

- 

- 

33.0 

32.4 

- 

- 

(G) 

28.6 

48.5 

47.4 

14.9 

- 

58.0 

54.3 

64.0 

55.7 

(R) 

48.7 

60.0 

59.2 

53.6 

« - 

74.2 

74.4 

79.7 

79.3 

(IR) 

- 

- 

76.7 

78.5 

- 

- 

- 

62.8 

74.3 

Method: 

A 

B 

A 

B 

A B 

A 

B 

A 

B 


1:7620 1:30480 


TRI-X 

459-3 


22.4 45.0 


11.6 31.3 


67.5 69.4 


A B 


TABLE 2. LINEAR REGRESSION ANALYSIS SUKNARY 


S kg/ha 


Film Type: COLOR COLOR IR TRI-X 

Frar.e: 457-25 456-25 459-25 

Plots (Bands) 


COLOR COLOR IR TRI-X 

457-3 455-3 459-3 


ALL(4,5,6) 

133.0 

110.0 

168.0 

149.0 

138.0 

149.0 

155.0 

130.0 

134.0 

112.0 

203.0 

172.0 

(G,R) 

163.0 

147.0 

164.0 

149.0 

- 


158.0 

126.0 

145. 0 

122.0 

- 


(B) 

212.0 

224.0 



- 

- 

123.0 

134.0 

- 

- 

- 

- 

(G) 

177.0 

153.0 

169.0 

170.0 

- 

- 

162.0 

144.0 

155.0 

125.0 

- 

- 

(R) 

166.0 

143.0 

162.0 

144.0 

- 

- 

160.0 

125.0 

149.0 

125. G 

- 

- 

(IR) 

- 

- 

218.0 

168.0 

- 

- 


- 

157.0 

131.0 

- 

- 

S0.(4,5,6) 

133.0 

88.9 

168.0 

133.0 

127.0 

91.1 

166.0 

162.0 

149.0 

£3.0 

233.0 2C3.0 

^(G,R) 

164.0 

163.0 

166.0 

139.0 

- 

- 

163.0 

145.0 

154.0 

141.0 


- 

(B) 

214.0 

247.0 

- 

- 


- 

103.0 

173.0 

- 

- 

- 

- 

(G) 

177.0 

155.0 

170.0 

128.0 

- 


169.0 

143.0 

153.0 

133.0 

- 


(R) 

162.0 

149.0 

163.0 

139.0 

- 


. 159.0 

142.0 

154.C 

141 .0 

• 

- 

(IR) 

- 

- 

223.0 

197.0 

- 

- 

- 

- 

161.0 

151.0 

- 

- 

S02(4,5,6) 

139.0 

134.0 

130.0 

95.4 

163.0 195.0 

135.0 

145.0 

121.0 

123.0 

132.0 

132.0 

& (G,R) 

159.0 

164.0 

155.0 

173.0 

- 

- 

125.0 

132.0 

110.0 

115.0 

• 

- 

A.R.(B) 

216.0 

218.0 

- 

- 

- 

- 

189.0 

195.0 

- 


- 

- 

(G) 

195.0 

171.0 

168.0 

220.0 

- 

- 

150.0 

151.0 

139.0 

140.0 

- 

- 

(R) 

165.0 

151.0 

148.0 

162.0 

- 

- 

117.0 

121.0 


1C3.0 

- 

- 

(IR) 

•” 


112.0 

111.0 

— 

— 

— 


141.0 

121.0 

— 

— 

Method: 
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B 

A 
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A 

e 

A 

B 

A 
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T/UJLE 3 


STANDARD DEVIATION SL’, ‘DIARY 


YioUl Cl.r. si f j j .Uion 

A porolleU'piped cliv.i.1 f Ication aloori llim wos written wliich, for eoth 
pixel, firtit eulcuKited yield from a specified repression equation and then 
placed that pixel In a yield class based on selected ranges. Pixels whose 
calculated yield did not fall in any of the specified ranges vvere designated 
unclassi fled, 

A subset of Inwge 62-fi56*"3 v/as classified using this method. Tlie 
regression equation of: 

Y - 1059. - 715000Xj 629000X2 - 7S3^i7Xj 

where: 

Y » predicted yield (kg/ho) 

X| « relative exposure of the green-sensitive layer 

X2 “ relative exposure of the red-sens I tivc layer 

X, *= relative exposure of the infrared-sensitive layer 

was used to calculate predicted yield for each pixel. This Image and regres- 

2 

slon equation combination was chosen to maximize R and minimize scale. The 
regression was derived from SO2 plots only (ILC, IL, 2LC, 2L, MC, M, HC, and H) 
and utilized 3 predictors (band band 5, and band 6). Additionally, 
plots for a particular treatment were merged prior to the regression (method 
B) . The sensor used was color Infrared film at a scale of 1:30,675. R and 
S may be found by reference to table 2 and table 3. 

Figure 1 graphically depicts the way In which the class intervals 
were chosen. The ranges In measured yields were plotted for each treatment 
type, and then, the ranges for the class Intervals were selected to correspond 
to the differences between the measured yields. 

Figures 2, 3» 5, and 6 are colored CRT displays of the classification 

and successive generalization routines. Figure 2 serves to illustrate the 
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graphics ovtM’lay proeetlure ami th« location of the plots In tlu* elassHled 
scene. Correspo'.dlng line prltucr maps may be found In the appendix. The 
classified area represents a scene Wil rows by 116 columns or ?16 M In the 
N-S dlrectloti by 178 H In the E-W diieetlon (each pixel Is approx Irma to ly 
1.53 M on 0 side.) for orientation purposes, the dirt read east of the plots 
Is depleted on the extreme right of the figures by a series of unelacsifled 
pixels In a vertical string. The class Intervals and their associated colors 
appear In the legend lioxes at the extreme right. Again, the yield ranges 
correspond to tlioso delineated in figure 1. 




tIElD K#/H«(XI00J 


JVMBOl 



KANOt 


*0 


% 


Figure 1. Measured Yield and Class Interval Ranges. 

9 


AnalytU of figure 3 indicates that, on a per pixel basis, there is 
considerable variation In yield even within a plot. Generally, higher yields 
•re associated with the controls, and lower yields arc associated with the 
SO2 damaged plots. However, more pixels were assigned to the ranges 600*1700 
and 3000*i«000 kg/ha than any other classes. These ranges seem to be probably 
outside the expected yields. Therefore, one must conclude that the rclotion* 
ship between yield and exposure is valid for a relatively .larroiv range. Per- 
haps a better equation could be found if a greater range of control plots 
could be found. 

figures 3, and S are the result of general izat ion routines. Pixels 
were reassigned to classes depending on the classes of their nearest neighbors. 
Figure seems to be the best compromise between Interprctabi I i ty and 
accuracy. However, this is a qualitative Judgement as no other plots were 
available for an independent measure of accuracy. 



Vv • r* 


I u 














CONCLUSIONS AND RCCO^^MCNDAT IONS 


The principal conclusion from this study must be that there is a 
definite relationship between film exposure ond yield in SO2 and acid rain- 
damaged soybeans. However, multiple linear regressions provide only a first 
approximation when a field situation is presented. Different stress agents 
seem to affect the relationship bct\ n exposure and yield in analogous but 
slightly different ways. Thus, quantitative description of the relationship 
requires extremely careful experimental design and meticulous attention to 
possible interactions between variables. While this study did not completely 
define the relationship betv/een exposure and yield, we feel that it is an 
Important first step in defining this relationship. 

It must also be concluded that as flight height increases, color 
Infrared film becomes the most useful sensor. This is most probably due to 
atmospheric interactions occurring in the blue region of the spectrum {.k - 
.5 pm). Black and white (Tri-X) film is of the least use at higher flight 
altitudes but is reasonably usable at lower ones. 

As one might expect, merging of plot data Improves the results, as 
does eliminating as many shadow and soil data points as is possible. 

Exposures from the red (chlorophyll absorption) band are most highly 

correlated with yield, with the green band producing redundant data. The 

Infrared and blue bands were found to be relatively small contributors at 

2 

high altitudes, but increased R In the low altitude data. 

From the above conclusions we recommend that if further study is to 
be made of this subject: 1 ) more plots be harvested such that the complete 

range of field conditions could be approximated, 2 ) each plot sho'ild be 
marked by a white panel, and survey data should be available for the location 
of each panel, 3 ) for efficiency, color infrared film should be flown at a 




16 


scale near 1:30,000, nullibond pholoyrapliy should be Investigated, J>) the 
multidate aspect of this study should be invest iyated, beginning with a fligiit 
at the end of July and continuing tlirough early September, and 6) it Is 
suggested that the overall classification accuracy would be Improved If a 
maximum likelihood classification was performed on the data prior to the 
yield classifier such that only soybean pixels would be classified using the 
yield class! f ier. 
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